IC/97/113 
hep-ph/9708403 



Vacuum oscillations of solar neutrinos: correlation between 
spectrum distortion and seasonal variations 



S. P. Mikheyev*'* and A. Yu. Smirnov^'* 



Oh' 
(D 

^ ' * Laboratory Nazionale del Gran Sasso dell'INFN 

o\ ; 

(N . 1-67010 Assergi (L'Aquilla), Italy 

^ I International Center for Theoretical Physics 

O ■ P- O. Box 586, 34100 Trieste, Italy 

00 

o 

ix. 

Abstract 

X 

' Long length vacuum oscillations solution of the solar neutrino problem is dis- 
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cussed. We show that there is a strict correlation between a distortion of the 
neutrino energy spectrum and an amplitude of seasonal variations of the neutrino 
flux. The slope parameter which characterizes a distortion of the recoil electron en- 
ergy spectrum in the Super-Kamiokande experiment and the seasonal asymmetry 
of the signal have been calculated in a wide range of oscillation parameters. The 
correlation of the slope and asymmetry gives crucial criteria for identification or 
exclusion of this solution. For the positive slope indicated by preliminary Super- 
Kamiokande data we predict (40 - 60) % enhancement of the seasonal variations. 



1. Long length vacuum oscillations of neutrinos on the way from the Sun to the Earth are 
considered as a viable solution of the solar neutrino problem [|], ||, ^ || (see P, 0, ||, || for 
latest analysis) There are two key signatures of this solution: (i) distortion of the neutrino 
energy spectrum and (ii) seasonal variations of the fluxes 0. 

(i) . The distortion follows from a dependence of the oscillation survival probability, P, on 
the neutrino energy E: 

P = l-sm'2esm'(^^^^ . (1) 

(We consider mixing of two neutrinos.) Here Am^ = — ml is the neutrino mass squared 
difference, 6 is the mixing angle and L is the distance between the Sun and the Earth. A 
variety of distortions of the boron neutrino spectrum is expected depending on values of Am^ 
and sin^ 29. 

(ii) . Seasonal variations are stipulated by ellipticity of the Earth's orbit. The flux of 
neutrinos at the Earth, F, can be written as 
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F = Fo[-^) -PmE), (2) 

where Fq is the flux at the astronomical unit Lq. The distance between the Sun and the Earth 
at a given moment t equals: 



L^Lo(l + ecos^j , (3) 



e = 0.0167 is the eccentricity, T = 1 year. The variations are expected both due to the 
geometrical factor, and due to change of the oscillation probability |5|. Depending on 

values of the oscillation parameters Am^ and sin^ 29 one may get an enhancement or damping 



of the geometrical effect |]TT], |12|, |13| or even more complicated variations of signals. 

In this Letter we point out that there is a strict correlation between a distortion of the 
boron neutrino energy spectrum and seasonal variations of the boron neutrino flux. This gives 
crucial criteria for identification or discrimination of vacuum oscillations solution. 



2. The correlation between the distortion and time variations can be immediately seen 
from the expression for the oscillation probability. Indeed, Eq. (|l|) gives 

dP _ dP E _ PE 

dL^'dE'L ~ ~' 

^ Large mixing angles implied by this solution are however disfavored by the data from the SN87A |]l0||. 
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where 

sAE) = (5) 

is the slope parameter which characterizes the distortion of the neutrino spectrum at a given 
energy E. If spectrum distortion has a positive slope, > 0, then according to Eq. (§) one 
gets dP/dL < 0, i.e. with increase of the distance the survival probability decreases. The 
change of the flux is in the same direction as due to the geometrical factor. Therefore, for 
a positive slope the vacuum oscillations enhance seasonal variations. For a negative slope, 
Sy < 0, Eq. (H) gives dP/dL > 0, that is, the probability increases with distance. In this case 
the oscillations weaken seasonal variations due to pure geometrical factor. 

Let us find the correlation explicitly. From Eqs. (|l],0) we get the change of neutrino flux 
with distance: 

dF F ( dP E\ 

Here the first term is due to the geometrical effect and the second one corresponds to a 
change of the probability with distance. For a positive slope both terms have the same sign 
in accordance with previous discussion. Introducing the seasonal asymmetry: 

A. = 2^, (7) 

where AF is the difference of the averaged fluxes during the winter and the summer and F is 
the averaged flux during the year, we can define the quantity 

A„ - AO 



with A'l = Ay{P = 1) being the asymmetry without oscillations. The asjTiimetry gives the 

variations due to the oscillations in the units of pure geometrical effect. From we get final 

relation between the slope parameter and the oscillation asymmetry r^: 

2 

Su = —Ty. (9) 

Let us stress that this relation does not depend on oscillation parameters Am^, sin^ 26 at least 
in the lowest order on the eccentricity e. 

3. The correlation described by (|^, |^) holds for fixed neutrino energy. In real experiments 
the energy spectrum of the recoil (produced) electrons is measured, and moreover, the inte- 
gration over certain energy intervals of neutrinos as well as electrons (due to finite energy 
resolution) takes place. This modifies the correlations. 
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Let us consider a manifestation of the correlation in the Super-Kamiokande experiment 
T^ , |T^ . For a majority of relevant values of Am^ and sin^ 26 the whole observable part of the 



boron neutrino spectrum (6 - 14 MeV) is in regions with definite sign of dP/dE. Integrations 
over the neutrino energy and on the electron energy weighted by the energy resolution function 
lead to a strong smoothing of the observable distortion of the recoil electron spectrum. The 



distortion can be well characterized by a sole slope parameter Se [16| defined as 



^^i?0 + SeTe, (10) 

where Nose and Nq are the numbers of events with and without oscillations correspondingly, 
-Ro is a constant, Tg is the recoil electron energy in MeV, Se is in the units MeV~^. There is 
some change of the spectrum distortion during the year due to variations of distance L |TT . 
However, this effect is small and can be neglected in the first approximation. 
The slope parameter is determined in the following way: 

(i) We calculate expected numbers of events, N^^^, for different values sin^ 29 and Am^ in 
the energy bins with ATg = 0.5 MeV from 6.5 MeV to 15 MeV (as in the Super-Kamiokande 
presentation of the data): 

KsciL) = j^^'^^^ej dT'JiTlX) l^dE^^^^FiE)PiE,L,Am\e), (11) 

where /(Tg, T^) is the energy resolution function, da{E, T'^/dT^ is the differential cross-section 
of the ve— scattering. 

(ii) Similar numbers of events, A^q, have been calculated in absence of oscillations: A'^q = 
A^*(P = 1). Then ratios N^^^/Nq have been found for each bin. 

(iii) The fit of the histogram NI^^/Nq by the function (|10D gives Sg (fig. 1 a). 



We will characterize the seasonal variations by the summer-winter asymmetry: 

Here ? ^sp^ ^ a are the numbers of events detected from November 20 to February 

19, from May 22 to August 20, February 20 to May 21, from August 21 to November 19 
respectively: 

Ny=fdtY.NUm) iV^W,S,SP,A). (13) 
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Similar asymmetry is calculated in absence of oscillations, i. e due to pure geometrical factor: 
Al = A,{P = 1). 

Finally, we introduce (as in (^) the signal asymmetry 

A^ 



'-e 



e 



which characterizes the asymmetry due to oscillations in the units of geometrical asymmetry. 
Positive Tg corresponds to an enhancement of the geometrical effect. 

4. The correlation between the spectrum distortion and seasonal variations appears as the 
correlation between the slope parameter Se and the signal asymmetry (fig. 1-4). 

In fig.l we show the dependence of the slope and asymmetry on Am? for different values 
of sin^ 26. Both the slope and the asymmetry increase with mixing angle, whereas positions 
of maxima and zeros (in Am^ scale) do not depend on 9. Note that zeros of the asymmetry 
are shifted with respect to zeros of the slope by approximately 8% towards bigger values of 
Am? . Maximal positive asymmetry is also shifted towards bigger Am^ (about 25% for the 
first maximum), however this difference diminishes for next maxima. In contrast, first negative 
maximum of asymmetry is shifted towards smaller Am^ . These features are related to the 
integration over the neutrino and electron energies. Maximal asymmetry due to oscillations 
can exceed the geometrical asymmetry (|re| > 1). 

Fig. 2 shows the slope - asymmetry plot. The points correspond to different values of 
Am^ and sin^ 26. Scattering of the points is related to the integration over the interval of the 
neutrino energies. According to fig. 2 the correlation can be approximated by 

Se = (0.03 - 0.10) MeV"^ . (15) 

Note that for the neutrino spectrum the slope is larger: according to (|^) Si, = 0.2 MeV^^rj, 
(r^ ~ re) for E = 10 MeV. 

Iso-asymmetry (a) and iso-slope (b) lines in Am^, sin^ 2^^ plot are shown in fig. 3. The 



iso-asymmetry plot (fig. 3a) is similar to the iso-plot of the near-far asymmetry found in |12 



Preliminary Super-Kamiokande data for 306 days [|r^ show a positive slope: Se ~ (0.013 — 
0.017) MeV~^. According to figs. 1, 3, in the region of the best fit of the integral data 
(Am^ = (5 — 8) ■ 10~^^ eV^) one expects for these values of the slope the positive asymmetry 
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Te = 0.4 — 0.6 . That is, an enhancement of the seasonal variations by (40 - 60) % should be 
observed. 

In fig. 4 we show real time seasonal variations of the fiux for different values of the 
neutrino parameters. Dashed-dotted line corresponds to very weak variations due to oscillation 
effect and therefore practically coincides with geometrical effect. Two other lines correspond 
to enhancement and suppression of the geometrical effect by oscillations. Solid line shows 
variations expected from the observed slope. Also shown are the preliminary results from the 



Super-Kamiokande (306 days) experiment |15| . Obviously, present statistics does not allow 



one to make any conclusion. Significant result can be obtained after 4-5 years of the detector 
operation. 

In larger interval of the recoil electron energies (from 5 MeV to 15 MeV) which will be 
accessible soon a description of the distortions by only one parameter may not be precise. 
Moreover, seasonal variations differ at different energies. In this case one can divide the 
interval into two parts: e.g., (5 - 8) MeV and (8 - 15) MeV and study the correlations in these 
intervals separately. A comparison of the effects will give the cross-check of the results. 



Similar correlation picture can be obtained for SNO experiment [[I^] where the effect is 
expected to be even more profound due to weaker averaging. 

Acknowledgment : 

Authors are grateful to E. Kh. Akhmedov and E. Kearns for discussion. 

References 



[1] 
[2] 
[3] 
[4] 
[5] 
[6] 
[7] 



B. Pontecorvo, Zh. Exp. Theor. Fiz. 53 (1967) 1717. 

V. Gribov and B. Pontecorvo, Phys. Lett., 28B (1969) 493. 

I. Pomeranchuk, in 0. 

V. Barger, R. J. N. Phillips and K. Whisnant, Phys. Rev. D24 (1981) 53^ 
S. L. Glashow and L. M. Krauss, Phys. Lett. B190 (1987) 199. 
P. I. Krastev and S. T. Petcov, Phys. Rev. D53 (1996) 1665. 
E. Calabresu et al., Astropart. Phys. 4 (1995) 159. 



6 



[8] Z. G. Berezhiani and A. Rossi, Phys. Rev. D51 (1995) 5229; Phys. Lett. B367 (1996) 
219. 



[9] N. Hata, P. Langacker, piep-ph/ 9705339] , (to be published in Phys. Rev. D) 



[10] A. Yu. Smirnov, D. Spergel, J. N. Bahcall, Phys. Rev. D49 (1994) 1389. 
[11] P. I. Krastev and S. T. Petcov, Nucl. Phys. B449 (1995) 605. 
[12] B. Paid et al., Phys. Rev. D55 (1997) 1353. 
[13] G. L. Fogh, E. Lisi and D. Montanino, BARI-TH/276-97, |hep-ph/9706228 



[14] Y. Totsuka, talk given at the Int. Symposium on Lepton and Photon Interactions (Ham- 
burg, July 1997). 

[15] K. Inoue, Talk given at the Int. workshop TAUP-97, Gran-Sasso, Italy, September 7-11 
(1997). 

[16] W. Kwong, S. P. Rosen, Phys. Rev. D51 (1995) 6159; A. Yu. Smirnov, Talk given at 17th 
International Conference on Neutrino Physics and Astrophysics (Neutrino 96), Helsinki, 
Finland, 13-20 Jun 1996, |hep-ph/961143^ . 



[17] G. T. Evan et al., Sudbury Neutrino Observatory Proposal, Report SNO-87-12 (unpub- 
lished) . 



7 




J I I I I I L 



J I I I L 



10 



-11 



10 



-10 



10 

Am^ (eV^) 



Figure 1: The dependence of (a) the slope parameter Sg (in MeV and (b) the signal 
asymmetry rg on Am^ for different values of sin^ 29. Solid , dashed, dash-dotted, and dotted 
lines correspond to sin^ 29 — 1.0, 0.75, 0.5, and 0.25 respectively. 
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Figure 2: The slope - asymmetry plot. The points (se — r^) correspond to different values of 
Am^ (200 points between 10"" and 10~^ eV^), and sin^ 26 (40 points between 0.025 and 1.0). 
Calculations have been done on grid 200 x 40. 
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Figure 3: a). Iso-asymmctry lines in Am^, sin^ 2^^ plot. Solid lines correspond to positive 
asymmetries: = 0.19, 0.38, 0.57, 0.76, 0.95 (from outer to inner lines). Dashed lines 
correspond to negative asymmetries: —Ve = 0.28, 0.56, 0.84, 1.12, 1.4. 

b). Iso-slope lines in Am^, sin^ 29 plot. Solid lines correspond to positive values of the slope 
parameter: Sg — 0.011, 0.022, 0.033, 0.044, 0.055 (from the outer to inner lines). Dashed lines 
correspond to negative slopes: -s^ = 0.013, 0.026, 0.039, 0.052, 0.065. 
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Figure 4: Real time seasonal variations of the flux for different values of the neutrino parame- 
ters: Am^ = 5.82- 10-11 eV^ sin^ 2^ = 0.90 (solid line); Am^ = 1.06-10-1° eV^, sin^ 2^ = 0.95 
(dashed line); Am^ = 7.85 ■ lO-^eV^ sin^26' = 0.75 (dash-dotted line). Dashed-dotted line 
coincides practically with the one due to pure geometrical effect. Also shown are preliminary 
results from the Super-Kamiokande (306 days) experiment pi|. 
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